Shot noise in the self-dual Interacting Resonant Level Model 
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By using two independent and complementary approaches, we compute exactly the shot noise 
in an out-of-equilibrium interacting impurity model, the Interacting Resonant Level model at its 
self-dual point. An analytical approach based on the Thermodynamical Bethe Ansatz allows to 
obtain the density matrix in the presence of a bias voltage, which in turn allows for the computation 
of any observable. A time-dependent Density Matrix Renormalization Group technique, that has 
proven to yield the correct result for a free model (the Resonant Level Model) is shown to be in 
perfect agreement with the former method. 



The study of transport in out of equilibrium, strongly 
interacting systems, is of crucial importance for techno- 
logical applications. It also addresses some of the most 
challenging fundamental questions, from the possibility 
of fluctuations theorems out of equilibrium to the 
time evolution of quantum many body entanglement . 

While experimental progress in this area has been swift 
and steady (see e.g. [sl-Q), the theory has been held 
back by considerable technical difficulties. The physics 
of interest occurs usually in non perturbative regimes, 
where analytical methods are few and limited. Numerical 
approaches require real time simulations, which up to 
recently had been notoriously difficult. 

Yet the situation is changing, especially in the case of 
quantum impurity problems. Extensions of the Bethe 
ansatz to study transport properties have been proposed 
[1,13, while time dependent DMRG (td-DMRG) calcula- 
tions have become quite reliable. This led to the break- 
through in [lo| where an exact, extremely non linear I-V 
characteristic was obtained analytically for the interact- 
ing resonant level model (IRLM) in the scaling limit, and 
recovered with remarkable accuracy in td-DMRG simu- 
lations on a lattice model. The results can now be used 
as a benchmark for approximate approaches. 

Interesting as it may be, an I-V characteristic does not 
capture the whole physics of a problem. In fact, the theo- 
retical and experimental activity has concentrated lately 
on fluctuations, embodied at vanishing temperature by 
the shot noise. Current fluctuations can be related to 
the charge of the "elementary" quasiparticles, and their 
study is akin to a subtle "many-body spectroscopy" . 

While the free case has been studied in details 0, fill . 
[T2j . progress on the interacting case has been very sparse. 
Tackling the shot noise certainly increases the level of dif- 
ficulty. On the analytical side, while the community has 
slowly accepted the idea that I-V characteristics could 



be obtained using the Bethe ansatz (thanks in part to al- 
ternate derivations, using all orders perturbation theory 
[l3| , or mappings onto an equilibrium system [l^ ) , there 
is, to our knowledge, only one published exact shot noise 
calculation [l^, which has not yet been confirmed. On 
the numerical side, extracting the noise from td-DMRG 
requires looking into the power spectrum of two point 
correlation functions, a quantity found to be strongly af- 
fected by finite size effects. 

Nevertheless, it has been realized very recently that the 
shot noise could indeed be reliably extracted from real 
time simulations, the finite size effects being controllable 
partly by analytical arguments, and partly by linear finite 
size extrapolations [la]. We use this idea to study the 
noise in the IRLM at the self dual point. Extending and 
adapting the ideas in [l^ we are able to obtain this noise 
analytically in the scaling limit, while we are also able 
to carry out the proposal of [l6j to determine this noise 
numerically. These two are in excellent agreement. 

The IRLM describes a single fermionic level (d^) that is 
coupled to two leads of spinless electrons playing the role 
of thermodynamical baths (reservoirs of charge and en- 
ergy) . The coupling involves a tunneling term that allows 
electrons to hop from the leads to the impurity, and an 
impurity/lead capacitive coupling (coulombic repulsion 
U). After the usual procedure of expanding on angular 
modes, unfolding the leads, and linearizing around the 
Fermi points, the Hamiltonian is given by: 

dxtPld^i'^ (1) 

Hb = (7iV^I(0)+ 72^1(0)) rf + H.c. (2) 
+ U{: V4V'i:(0)+ : 4iJ2-iO)) {dU - ^) + dU, 

where the baths are described by two ID right-moving 
fields ipa{x). The impurity /leads coupling involves the 
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tunneling amplitudes that we parametrize as: 

71 + i72 = 7V2 e'^/^ (3) 

In the following, we consider the (particle-hole symmet- 
ric) resonant case with impurity onsite energy ed = 0. 

The IRLM bears a duality symmetry exchanging large 
and small C/'s [l^. For an intermediate value of U (of 
order unity), it is self-dual, and enjoys an additional, hid- 
den, SU(2) symmetry that mix the two wires, with gen- 
erators J = h^l<J„K^h (cr's are Pauli matrices) ; as was 
shown in RefflOl. as a consequence, the out-of-equilibrium 
self-dual IRLM (sd-IRLM) bears a description in terms 
of dressed quasiparticles (qp's), that are the many-body 
modes diagonalizing the scattering on the impurity. 

Let us now sketch the main lines of the derivation. As 
usual in a scattering approach, one starts by identifying 
two classes of asymptotic states, incoming and outgoing, 
that correspond to states coming from the far left to- 
wards the impurity, and escaping to the far right thereof, 
respectively. Those states span Hilbert spaces and 
J^Qut- The effect of the impurity amounts to a linear 

map — > y^out that encodes the fate of an asymp- 
totic in state prepared in the far past when time-evolved 
to far future. As soon as there is an interaction, this 
linear map becomes a complicated many-body object in 
the electronic basis. Integrability of the IRLM en- 
sures that one can identify a basis for and J^out in 
terms of pseudo-Fock states built out of a finite number 

of quasiparticle modes (0), where is a rapidity 

parametrizing momentum, p = ^e®. Those modes di- 
agonalize the map TZ, in the sense that they cross the 
impurity without qp production, and that this property 
extends to any many-qp state [1^ . The only effect of the 
impurity is to change the qp index a. Formally, one has: 

Afi9) ^ R^,i9)A$-H9) (4) 

with R a scattering matrix. For the sd-IRLM, such a 
basis can be obtained via bosonization and a mapping to 
the anisotropic Kondo model (see [l^l for details). The 
total charge degree of freedom decouples from the prob- 
lem, and one is left with a single degree of freedom, the 
charge imbalance between the two wires. The qp's con- 
sist of a soliton and an antisoliton A± , and two breathers 
Aq and Ai. Importantly, those qp's fall into representa- 
tions of the aforementioned SU(2) symmetry: Ai} 
transform as a spin one, while Aq is a singlet. The charge 
imbalance Q = J dx[iplilji — iplip2) —2 J dxJ^ acts diago- 
nally on the modes: Q ■ Aa{0) = qaAa{9), with q±=zL2e 
and (70.1=0. Introducing the operator A^^ that destroys 
the qp Aa, the charge Q bears a simple representation 
in terms of the modes: Q = J2a I '^^^a^a{d)Al^{9). 
To complete the description, non-vanishing elements of 
the i?-matrix are R±± = Q-V cos^ F, R±^ = V sin^ F, 
Ri±=R±i=V^^^^ , i?ii = Q+7'cos(2F), and i?oo, with 
= n^^o.ii e-e,t^ . and Q{9) = -ie3(«-«B)p((^). 

e 3 -f 1 

The lead/impurity coupling results in the appearance of 
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FIG. 1: Finite size error of noise. The blue lines represent 
the analytical result Eq. (|6]) with symmetric coupling {T=^). 
The numerical data have been obtained for systems with cou- 
pling J' = 0.4 J and density-density interaction (7 = 2.0 J using 
td-DMRG. The system size varies from M = 48 to M = 72 
lattice sites. The difference of numerical and analytical data 
in the low voltage regime is proportional to the squared dif- 
ferential conductance and scales linearly with the inverse 
system size 1/M. 



an energy scale Tb = ^ ^ marking the crossover be- 
tween weak and strong hybridization regimes. 

Forcing the sd-IRLM out of equilibrium is achieved by 
imposing different chemical potentials /ii(2) = on 
incoming electronic states in wires 1 and 2, i.e. by cou- 
pling the system to ^Q^^ = \ dx{iplip^ - ipl'ip2)- 
In the qp basis, this operator acts diagonally on the 
modes A™, and the voltage Vsd just amounts to a doping 
of in modes. At zero-temperature and positive voltage, 
the groundstate is obtained by filling antisoliton states 
up to a voltage-dependent rapidity A = ln(2pF /m) where 
Pf oc Vsd is akin to a "Fermi momentum" (because of in- 
teractions, the proportionality constant is not one), with 
distribution function (9) that is determined by doing 
the thermodynamics for the gas of incoming antisolitons. 

The current operator counts the charge imbalance be- 
tween in and out modes, and reads (the factor i comes 

from simple charge counting): I = ^{Q^'^ — Q°^^), that 
can be rewritten solely in terms of in modes using (j4]): 

/ = J d9A'^i9)U^p{9)Al'-{9) (5) 

with n = i {Q — R*QR) and Qap = qaSafs- Averaging 
© yields /=(/)= 2 sin^ F J^^ d6'p^^° {0)T{9), with T = 
IPp. The zero frequency noise in the steady state, So = 
^/•^) — is obtained by averaging the square of 

The terms (Ajf (6')A|,'°(6')Ajf, (6'')A^;°(6'')) that appear 
have a simple expression at zero temperature Q , and the 
noise reads: 5*0 = 2 sin^ F d9 fVf' {9) [(l-^sin^ T)T{9) - 
2 sin^ F(T(0))^] . Simple algebraic manipulations using 
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FIG. 2: Linear extrapolation of the numerical shot noise data. 
The linear scaling of the low voltage hnite size error is ex- 
ploited to perform a linear extrapolation 1/M — >■ 0. We find 
nice agreement of numerical and analytical results. 

the scaling form [l^ of the current, / = Vsd/(-^), yield: 

5*0 = cos^r /+ isin2r(/- GV'sd) (6) 
with G = -M— the differential conductance. For self- 

OVSD 

containedncss, the scaling function defining the current 
is: /(x) = ^/(f), a^ 'f^\%/'l , Go = G(ysD = 0) = 

sin^ r 1 f( \ _ \p (-l)"(4n)! gn ^ „ VS 

h "^^^^ y-^) ~ l^n>0 n!r(3n+3/2)'^ iOi 0, A — ^2/3 

while fix) = Y.n>o ^^irW53W^^^ ^""^ 

In the high voltage limit, the Fano factor So /el ~ 
^+'^°^ ^ = ^j-^'^l.o is identical to that of free electrons 
tunneling through a double barrier with tunneling rates 
proportional to 7^ [l^] . Thinking of voltage as a high en- 
ergy cut-off, this is consistent with the picture of charge 
carrier being independent electrons [2^. At low voltage, 
for a symmetric junction F = ^ , introducing the backscat- 
tered current Ibs = GqVsd — I, one obtains Sq ~ 2e/Bs, 
i.e. Poissonian noise for charge 2e particles, consistent 
with antisolitons being the charge carriers at low energy. 

A completely different approach is based on the nu- 
merical simulation of the time evolution of the system. 
We use the td-DMRG [2l| to obtain the time dependent 
current fluctuations 

sit,t') = {Ai{t)Ai{t'))^, Ai{t) = i{t)-{i{t))^, (7) 

/(t)-e'^*ic-'^', (•)* = (*| • 1*), (8) 
for a tight-binding lattice version of the IRLM 
H = Ho + Hb, (9) 

Ha = -j[ Y cL+iC,„+ Y cL+iC,„} +H.C., 

rn— — JI/l m— 1 

= ^ {u{n„~^){nd-^)-J'{cld+R.c.)'j+Vgnd, 

with the creation, annihilation, and density operators in 
the leads (c], Cj, n- = CjCj) and on the level {d'' , d, 



nd = d''d). The total number of sites is given by M = A4i^+ 
Afp^-l-l. The time dependent fluctuations then allow for 
the calculation of the noise power spectrum [23l by means 
of a Fourier transform S{uj) = iRe dte''^*S{t,t'), 
where we now restrict ourselves to the low frequency 
limit w = 0. Note that the discrete nature of the leads 
results in a finite bandwidth 4J with cosine dispersion. 
We create an initial state by taking the ground state of 
H+Vs-d{Ni,—N-r)/2, where we added a charge imbalance 
operator to the system H in order to create different fill- 
ings in the left and right lead. We now perform a time 
evolution with respect to H and search for a stationary 
regime close to the impurity, for details see [22j. 

Within this approach we have to address the following 
time scales. First, initial oscillations of the current die 
out on a time scale tg which is typi cally proportional to 
the inverse of the resonance width [22, 123 ■ Second, at 
transit time = M/v{ the reflections at the boundaries 
reach the impurity and we have to stop our simulations. 
Third, after choosing a starting point t' ^ tg we have to 
simulate until imax < ^r, where tmax^i' has to be large 
enough to allow for a faithful calculation of the noise 
power spectrum. A detailed analysis of the finite size 
correction in the non- interacting case is provided in (l6| . 

Here, we set U = 2.0J, corresponding to the sd-IRLM 
as discussed in [l^l, and the coupling to J' = 0.4 J, while 
we operate in the resonant tunneling regime = 0. The 
total number of lattice sites varies from A/ = 48 to 72 
lattice sites, with Mr = Ml + 1. Different other setups 
have been considered, including the effective enlargement 
of the system using damped boundary conditions, which 
will not be presented in this work. For the numerical 
simulation within the DMRG projection scheme we set 
an upper bound to the dimension of the Hilbert space for 
each DMRG block to iVcut = 4000 states. 

As a first result wc compare the numerical data for dif- 
ferent system sizes to the analytical result in Fig.[Tl where 
we show zero-frequency shot noise as well as the finite size 
error of the numerical data, rescaled by the system size. 
In the low frequency limit, strong finite size effects have 
to be expected, th at g et mostly pronounced for small val- 
ues of the voltage [lal • Since the rescaled finite size error 
happens to collapse on a single curve in the low voltage 
regime, the numerical data can be linearly extrapolated 
to infinite system size in order to obtain results for the 
thermodynamic limit. Also we verify an analytical esti- 
mate for the finite size error S'num. — i^anaiyt. esc G^/M 
with G the differential conductance that has been given 
for the free case in [l^ . The strong deviations in the high 
voltage regime from this relation may be traced back to 
different sources: the approximative td-DMRG scheme 
introduces a cutoff error that gets especially pronounced 
for values of the voltage of the order of the bandwidth. 
Furthermore, to keep the numerical simulation feasible, 
one has to resort so small systems introducing finite size 
effects beyond the linear scaling. 

Nevertheless, the numerical results shown in Fig. [2l 
where we obtained data for the low voltage regime using 
linear extrapolation, show very nice agreement with the 
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FIG. 3: Back scattering Fano factor as a function of the back 
scattered current. The numerical data points have been ob- 
tained using the numerical shot noise data divided by the 
analytical back scattered current. The finite size error of the 
numerical results for shot noise leads to a diverging Fano fac- 
tor. The situation improves for the linearly extrapolated data, 
while we find a nice agreement of the analytical result with 
the G^-corrected data. For comparison we show the Fano 
factor in the non-interacting case. 



analytical results given by Eq. ^ with symmetric cou- 
pling, r = -J . The back scattering Fano factor Fgs = S/ Ibs , 
can also be obtained from the numerical data, Fig . [3l 
where we use the analytical result for the current [lO| . It 
fits nicely with the analytical result for F^s as long as fi- 
nite size effects can be neglected. However, the finite size 
offset at /bs leads to a strongly diverging Fano factor. 
In contrast, F^s remains finite even for very small values 
of Ibs, when obtained from the linearly extrapolated shot 
noise data. The deviations from the analytical result at 
small Ibs can be traced back to small absolute errors that 



get blown up in the limit Ibs 0. The very nice agree- 
ment of analytical result and G^-corrected data, even in 
the regime of very small Ibs , indicates that increasing the 
system size and adding more data points to the extrapo- 
lation procedure should improve the extrapolated result. 

In summary we have provided two different methods, 
an analytical approach within the framework of the ther- 
modynamic Bethe ansatz, and time dependent DMRG 
simulations on the lattice to obtain the noise correlations 
in the IRLM. Both methods show excellent agreement 
and provide benchmark results for other methods. Most 
strikingly our results show a strong enhancement of the 
back scattered Fano factor due to interaction effects. 

On the conceptual side, we believe our result further 
establishes the reliability of the Bcthc ansatz approach 
to transport pioneered in |15| . One of the objections to 
this approach sometimes raised is that it relies on the 
theory of excitations over the vacuum, and thus deals 
with fundamental objects - quasiparticles - which are not 
simply related to the bare electrons. In the present case, 
these objects are solitons, of charge \q±\ = \2e\, which 
arc made of combinations of particle hole pairs mixing 
the two wires. The results in the low energy limit give 
these quasiparticles their physical reality: they are the 
objects that tunnel in a Poissonian way at low voltage, 
and the Fano factor is a direct measure of their charge. 
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